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NLO QCD corrections to H, W and Z production via vector boson fusion have recently been calculated in 
the form of flexible parton level Monte Carlo programs. This allows for the calculation of distributions and cross 
sections with cuts at NLO accuracy. Some features of the calculation, as well as results for the LHC, are reviewed. 



1. Introduction 

Vector-boson fusion (VBF) processes have 
emerged as a particularly interesting class of scat- 
tering events from which one hopes to gain in- 
sight into the dynamics of electroweak symme- 
try breaking. The most prominent example is 
Higgs boson production via VBF, which is shown 
in Fig. 1. This process has been studied inten- 
sively as a tool for Higgs boson discovery [1-3] and 
the measurement of Higgs boson couplings [4] in 
pp collisions at the CERN Large Hadron Collider 
(LHC). The two scattered quarks in a VBF pro- 
cess are usually visible as forward jets and greatly 
help to distinguish these Hjj events from back- 
grounds. 

Analogous to Higgs boson production via VBF, 
the production of Wjj and Zjj events via vector- 
boson fusion will proceed with sizable cross sec- 
tion at the LHC. These processes have been con- 
sidered previously at leading order for the study 
of rapidity gaps at hadron colliders [5-7], as a 
background to Higgs boson searches in VBF [2,3], 
or as a probe of anomalous triple-gauge-boson 
couplings [8], to name but a few examples. 

In order to match the achievable statistical pre- 
cision in such LHC studies, the inclusion of QCD 
corrections in the predicted VBF cross sections 
is required. While NLO QCD corrections to the 
Higgs boson total cross section have been known 
for over a decade [9], NLO corrections to distri- 
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butions have become available only recently, both 
for Higgs boson production [10,11] and for 
and Iv production in VBF [12]. These new calcu- 
lations use the subtraction method of Catani and 
Seymour [13] to construct flexible parton level 
Monte Carlo programs for the calculation of NLO 
corrected distributions. In this talk we describe 



the work of Rcfs. [10,12]. 
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Figure 1. Feynman graphs contributing to the 
VBF process qQ qQH at (a) tree level and (b) 
including virtual corrections to the upper quark 
line. 

2. Elements of the calculation 

Because of the color singlet exchange in the LO 
VBF diagrams, any 0{as) corrections, where glu- 
ons are simultaneously attached to both the up- 
per and the lower quark line in Fig. 1, vanish 
identically. Hence it is sufficient to consider ra- 
diative corrections to a single quark line at a time. 
For Higgs boson production, the virtual diagrams 
reduce to simple vertex corrections, like the one 
depicted in Fig. 1(b). For W and Z production 
also box graphs appear in corrections to diagrams 
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where the final state vector boson is emitted from 
one of the quark lines. 

The soft and coUinear singularity structure for 
real emission corrections only depends on the 
color of the external partons and, therefore, is 
universal for all VBF processes. Consider, for ex- 
ample, the emission of a gluon from the upper 
quark line in Fig. 1(a), with momenta given by 

qiPa) + QiPb) ^ 9{pi) + q{p2) + Q{P3) + B{P) .(1) 

Here B = H,W,Z denotes the produced bo- 
son and we write the trec-lcvcl amplitude for 
the emission process as Mf = Mf{pa,Pi,P2]q), 
where q = pi + P2 — Pa is the four momentum of 
the vector boson V which is attached to the up- 
per quark line and which has virtuality = —q^. 
The singularities of the 3-parton phase-space in- 
tegral of can be absorbed into a single 
counter term 
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where 



■M'^g{pa,P2', q) is the Born ampli- 
tude for the lowest order process 

q{pa) + Qipb) ^ q{p2) + Q{pz) + B{P) , (3) 
evaluated at the phase-space point 



Pa = xpa 



P2 = Pi + P2 - (1 - X)pa , (4) 



with 

X = 



Pi ■P2 



(Pl + P2) ■ Pa 
Pi - Pa 



P2-Pa 



{Pl+P2)-Pa {Pl+P2)-Pa' 
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This choice continuously interpolates between the 
singularities due to final-state soft gluons {pi — > 
corresponding to x — » 1 and z ^ 1), coUinear 
final-state partons {pi \ \p2 resulting in pi ■ p2 
or X — *■ 1) and gluon emission coUinear to the 
initial-state anti-quark {pi (1 — x)pa and 
2; ^ 1). The subtracted real emission amplitude 
squared, {Mfl"^ - |A1^|sing) leads to a finite phase- 
space integral of the real parton emission cross 
section, and these; integrals are evaluated numer- 
ically in D = 4 dimensions. 

The singular counter terms are integrated ana- 
lytically, in Z) = 4 — 2e dimensions, over the phase 



space of the coUinear and/or soft final-state par- 
ton. 1/e terms proportional to the splitting 
function disappear after factorization of the par- 
ton distribution functions. The remaining diver- 
gent terms in the integral of Eq. (2) yield the con- 
tribution (we are using the notation of Ref. [13] 
and Cf = 4/3) 
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The 1/e^ and 1/e soft and coUinear divergences 
cancel against the poles of the virtual corrections. 
For Higgs boson production they are depicted in 
Fig. 1(b) and are given by a simple vertex cor- 
rection only, which is proportional to the Born 
amplitude. For the more general case of W and 
Z production, the virtual amplitudes involve box 
contributions and have a much more complex 
structure. However, by isolating the 1/e^ and 1/e 
poles of the Passarino-Veltman functions, one can 
show that the divergent contribution to the sum 
of all virtual graphs is again proportional to the 
overall Born amplitude. 
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where Ady is finite. Thus, the interference be- 
tween the Born amplitude and the virtual correc- 
tion amplitude, 2K [A^\/A^|j], exactly cancels the 
divergent terms in Eq. (7). The remaining inte- 
grals, involving the finite remainder 25R[A^yA1^], 
are performed numerically in D = A dimensions. 

3. Predictions for the LHC 

The; calculations discussed above have been im- 
plemented in the form of a NLO Monte Carlo 
program. This permits the determination of ar- 
bitrary infrared and coUinear safe distributions 
and cross sections within cuts. In order to recon- 
struct jets from the final-state partons, the kr- 
algorithm [14] as described in Ref. [15] is used. 
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Figure 2. Total Higgs boson production cross sec- 
tion as a function of the Higgs boson mass, within 
the cuts described in the text at LO (dotted hne) 
and at NLO. The NLO curves are shown for se- 
lection of the tagging jets as the two most ener- 
getic jets (E method) or as the two jets of highest 
transverse momentum (px method). 

with resolution parameter D — 0.8. As a de- 
fault, electroweak parameters are determined in 
the scheme, with mz — 91.188 GeV, mw — 
80.419 GeV and the measured value of Gp as our 
electroweak input, from which we obtain uqed — 
1/132.51 and sm^ Ow = 0.2223, using LO elec- 
troweak relations. The decay widths are then cal- 
culated as Tw = 2.099 GeV and Tz = 2.510 GeV. 
We use the CTEQ6M parton distribution func- 
tions (PDFs) [16] with as{mz) = 0.118 for all 
NLO results and CTEQ6L1 parton distributions 
for all LO cross sections. For further details, see 
Refs. [10,12]. 

In the following we consider Hjj, Wjj and Zjj 
cross sections within generic cuts which are rel- 
evant for VBF studies at the LHC. We ask for 
events with at least two hard jets which are re- 
quired to have 

PT, >20GeV, |%|<4.5. (9) 

Here yj denotes the rapidity of the (massive) jet 



momentum which is reconstructed as the four- 
vector sum of massless partons of pseudorapid- 
ity |ry| < 5. The two reconstructed jets of 
highest transverse momentum are called "tagging 
jets" and are identified with the final-state quarks 
which are characteristic for VBF processes. 

We calculate cross sections for decays Z — > 
£^£^ and W ^ iv into a single generation of lep- 
tons. In order to ensure that the charged leptons 
are well observable, we impose the lepton cuts 



PTi > 20 GeV , 



<2.5, ARje>OA, (10) 



where Rji denotes the jet-lepton separation in the 
rapidity-azimuthal angle plane. In addition, the 
charged leptons are required to fall between the 
rapidities of the two tagging jets, 



Uj.Tnin ^ ^£ ^ yj.rnax ■ 



(11) 



For the case of Higgs boson production, the ap- 
proximately massless Higgs boson decay products 
T+T^, bb or 77 play the role of the leptons and 
we impose the above lepton cuts on them, except 
that ARji > 0.6 is used in the Higgs boson case. 
Since the branching ratio, B, for each of these fi- 
nal states depends on the respective coupling to 
the Higgs boson, we show Higgs boson cross sec- 
tions renormalized to a branching ratio of 100% 
for the selected final state, i.e. the cross section 
within cuts is multiplied by an overall factor 1/ B. 

Backgrounds to vector-boson fusion are signif- 
icantly suppressed by requiring a large rapidity 
separation of the two tagging jets, i.e. we require 



^Vjj = Ivji -yh \ > 4 



(12) 



The Higgs boson production cross section, 
within the above cuts, is shown in Fig. 2 as a 
function of the Higgs boson mass. QCD correc- 
tions increase the LO cross section (dotted line) 
slightly, by about 3 to 8%, depending on the Higgs 
boson mass. In addition to our default of identi- 
fying the tagging jets as the jets of highest trans- 
verse momentum (px method). Fig. 2 also shows 
the NLO result which is obtained when defining 
the tagging jets as the two most energetic jets in 
a given event {E method, dashed line). From a 
consideration of tagging jet pT-distributions, the 
Pt method appears to be somewhat more stable 
against scale variations, however [10]. 
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Figure 3. Scale dependence of the Z 
cross section in VBF at the LHC. 
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The typical scale dependence of the NLO VBF 
cross sections is demonstrated in Fig. 3, where 
the Zii cross section within the cuts of Eqs. (9- 
12) is shown as a function of the renormalization 
scale, [iji = ^mz, holding the factorization scale 
fixed at fj,p = mz (dashed curve), as a function 
of ^ip — £.niz for /i/j = mz (dot-dashed curve) 
or when varying both. In all cases one finds scale 
variations at below the 2% level (when allowing 
variations between ^ = 0.5 and ^ = 2) which 
is minute compared to the LO variation (dotted 
curve). This result holds for Higgs, W and Z 
production as well as for the scale choice fi = £,Q, 
where Q denotes the virtuality of the t-channel 
weak bosons in the various VBF processes. 

The NLO QCD corrections to observable VBF 
cross sections are modest, amounting to an in- 
crease of less than 10% for ^ = 1 in the ex- 
amples discussed above. Similarly modest QCD 
corrections are found for most distributions also. 
One example is shown in Fig. 4, where the tran- 
verse momentum of the hardest tagging jet in 
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Figure 4. Transverse-momentum distribution of 
the highest-pT jet in production at the LHC. 
In panel (a) the differential distribution is shown 
at LO and NLO for the scale choice hf — M_r = 
mw [M scheme). In panel (b), we show the ratios 
of the NLO differential cross section in the M 
scheme (solid black line), of the LO one in the 
M scheme (dashed red line) and of the LO one 
in the Q scheme (blue dotted line) to the NLO 
distribution for the scale choice hf — A*-R = Qi- 
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events is given. A compaxison of the LO 
(dashed line) and NLO distribution (sohd hne) in 
panel (a) exhibits modest shape changes. They 
are more clearly exposed in Fig. 4(b) where the 
ratio of distributions 



RipT) = 



(13) 



is shown, i.e. we compare to the default scale 
choice ^ ~ Q at NLO. Using the W mass as 
the renormalization and factorization scale in- 
stead introduces changes below the 2% level at 
NLO (solid line) which again points to the sta- 
bility of the NLO predictions. At LO the differ- 
ential cross section for this scale choice is sup- 
pressed by 10% and more at low Pyt^gi ^^"^ 
enhanced by up to 10% at high px, i.e. there is 
a substantial shape change from the LO to the 
NLO prediction. These changes are somewhat 
less pronounced when choosing fi = Q for the LO 
calculation. 

4. Conclusions 

The one-loop QCD corrections to Higgs, W 
and Z production via VBF in hadronic collisions 
are now available in the form of a flexible NLO 
parton-level Monte Carlo program. The NLO 
corrections are modest for cross sections as well 
as differential distributions at the LHC, rarely ex- 
ceeding the 10% level. In a variety of distribu- 
tions one finds shape changes at the 10% level. 
These changes need to be taken into account in 
precision studies of VBF processes at the LHC. 
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